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We show that a recent measurement of surface superconductivity in UPt3 (Keller et. al, Phys. 
Rev. Lett. 73, 2364 (1994)) can be understood if the superconducting pair wavefunction is sup- 
pressed anisotropically at a vacuum to superconductor interface. Further measurements of surface 
superconductivity can distinguish between the various phenomenological models of superconducting 
UPt 3 . 
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Superconductivity in hexagonal UPt 3 exhibits many properties unique among superconductors. The most striking 
of these is the pressure temperature magnetic field phase diagram that shows several superconducting phases . 
For example, a graph or the upper critical field for fields in the basal plane (H^) versus temperature exhibits a kink 
at a temperature T* §j. This kink can be explained in terms of a generalized Ginzburg-Landau model in which the 



order parameter has two complex components, 771 and r\i ]5|-|8|. For temperatures below T* , one of these components 

n 
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orders at Hf 3 , whereas for temperatures above T* the other component orders at if" 



Recently, a detailed study of surface superconductivity in UPt3 by Keller et. al. has produced puzzling features 
not observed in other superconductors |)|,[l0) (see also In particular, the ratio of the upper critical field for 

surface superconductivity (H C3 ) to that for bulk superconductivity has an unusual temperature dependence p0[ . 
Furthermore, different surfaces of the UPt3 whisker were observed to have different values of H C3 [|||l0|]. It is argued 
here that the reason for this unusual temperature dependence is that one of the components for the order parameter 
(say ?7i) is suppressed at the surface (and hence cannot sustain surface superconductivity); the transition at H C3 is 
thus to a superconducting 772 state, both below and above T*. Order parameter suppression at a superconductor to 
vacuum interface is known not to occur in conventional isotropic superconductors |13] but is expected on theoretical 



grounds (although not yet observed experimentally) for some unconventional superconductor symmetries [[13 14 1. We 
believe the observed behavior of H C3 /H C . 2 in UPt3 is the only evidence to date of an anisotropic suppression of the 
superconducting order parameter at a vacuum to superconductor interface. This interpretation of the measurements 

S of Keller et. al. imposes strong constraints on existing phenomenological models of superconductivity in UPt 3 
and indicates that such measurements can provide a new technique for determining the symmetry of anisotropic 
r 7^ • superconductors. 

In this article we will initially consider surface superconductivity for the _B-REP Ginzburg-Landau models proposed 
by by Machida and Ozaki || and by Hess et. al. || and improved by Sauls [|| (an important element of these models is 
the hexagonal symmetry breaking due to the basal plane antiferromagnetism [ [^[Lsl ). We then show how the inclusion 
of symmetry dependent surface contributions to the free energy is essential for understanding the experiments. Finally, 
we show how surface superconductivity studies can distinguish between these models, the ID-REP model [^8|, and 
the AB model (T|]. 
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For the odd parity Em and E\ u models the gap function is described by the pseudo-spin pairing gap matrix 
A(r, k) — i[rji(r)di(k) ■ a + 772(^)^2 (k) • <?]cy where Ui are the Pauli matrices. We will assume that di(k) — di(k)z; this 
form of the gap function is required in the E 2u model to account for the gap nodal structure implied by experiment 
H . For the even parity E 2g and E\ g models the form of the gap matrix is A(r, k) = ia y [^i{r)il)i{k) + rj2(r)ip2(k)]. 
The free energy density for these E-KEP models is 

F = a (T - T c )rf- rf - 7(^2 1 2 - M 2 ) + 4lAr/ 2 | 2 + ^\D lVl \ 2 

4\DM 2 + ^\D zm \ 2 + K 2 {D i r, i ){D j r 1j y + k 3 (A%)(A%)* + h 2 /8ir - h ■ H/An (1) 

where D = d — (i2e/Tic)A, h = d x A, H is the applied field, i and j refer to indices x (and 1) and y (and 2), 
and Kr = k/(1 ± Si). The convention a^rji = —771 for even parity representations and 03771 = 771 for the odd 
parity representations fixes the axes and is a reflection in a plane normal to the a crystallographic axis (these 
conventions allow the discussion of surface superconductivity to occur on the same footing for the even and odd parity 
representations). For the calculation of H C2 and H C3 higher order terms in rji and 772 in Eq. ^ can be omitted. The 
coefficients 7 and Si arise from the coupling to the antiferromagnetic moment. The gradient terms with coefficients 
K2 and K3 give rise to anisotropy in 20 1. This has been observed to be of the order 0.01if° 2 b |J. It has also 

been argued by Sauls || that these terms are small relative to K\ in the E^u model. We will therefore solve for H C3 
to first order in A = (K2 + k,^,)/2k\. The linearized Ginzburg-Landau equation given by the free energy in Eq. [l] is 
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->_ I ^1-23-°^ + K i D l + K 4 D l ~1 K 2 D x D y + n 3 D y D x »_ 

'"' 1 K3 D x D y + K2 D y D x Kl+23 Dl + K+D 2 x + K+Dl+ 1 ]l ' 



where Ki±23 — + k 2 + k 3 . 

We will consider the geometry in which the superconductor occupies the half space y > 0. The surface normal will 
be taken to be along the a (or y) direction and the antiferromagnetic moment along the a* (or x) direction. The field 
will be chosen to lie along the a or a* directions. This geometry will allow for a comparison to be made with the 
measurements of Keller et. al. [9 10 1. 



The off diagonal terms in Eq. |2| containing n 2 and k 3 vanish in the determination of the upper critical field for the 
field orientations chosen above when there is no surface present. In the presence of a surface normal along a it can 
be shown that to first order in A these terms have no effect. In this limit the Ginzburg-Landau equations decouple 
for the 771 and r) 2 components. When the field lies along a then the upper critical field (denoted H£ 2 ) is given by 



\Jk+4[H^] 2 = -a + 7 , m ^ m = 



«r«4 (l + A)[i? c a Ji = -a -7, 771^0 7/2 = 0. (3) 

where a = a^(T — T c ). If 7 > and 5\ + S 4 > or 7 < and 61 + 64 < then the two solutions will cross at a 
temperature T*. This crossing gives rise to the kink in H® 2 (T) that is observed in experiment [Q. For the field along 
the a* direction the upper critical field for surface superconductivity (denoted H^ a where a refers to the surface 
normal and a* to the field direction) can be found using the exact solution for conventional superconductors Q ; this 
yields 



+ K+(l + \)fi 2 [H2f} 2 = -a + 7, V2^0 771 = (4) 



«7 

V^/J 2 [.ff c a 3' a *]i = -a - 7 771 ^ 772 = (5) 

where /x 2 = 0.59010. An immediate consequence is that H^ a will have a kink where the two solutions cross. This 
kink will occur at the approximately T* (the difference is of order A). Note that H^ a /H® 2 — 1/fi 2 + 0(A) for all 
temperatures in this model. Keller et. al. PJlO|] have observed that this ratio is strongly temperature dependent. 
This model cannot account for the observed behavior. 

Qualitatively different behavior from that presented above arises when the surface free energy is considered |l4j . For 
conventional isotropic superconductors with order parameter ip, the surface free energy takes the form g J sur j: ace \ip\ 2 dS. 
For g > this term suppresses surface superconductivity. However, superconductor to insulator or vacuum boundaries 
of isotropic superconductors are well described by g — |H| . For anisotropic superconductors this is not necessarily 
the case. For the E'-REP models and the geometry defined above the surface free energy density has the form 

Fsurface = gi\m\ 2 + 92\m\ 2 - (6) 

In the presence of this surface free energy density the boundary conditions are 

- (n^Dyr/x + K 3 D x ri 2 ) = .91771 
-(ki+23-D»»?2 + n 2 D x rn) = g 2 r] 2 . (7) 

Assuming that the coefficients are field independent one can find H^ a by using the method similar to that used for 
an isotropic superconductor [l^] . For convenience we set A = in the remainder of this article (the corrections O(A) 
are easily obtained). The solution is given by 



where i refers to the solutions ffx = (1,0) or t/ 2 = (0,1), h = [H^J H^f]\ /2 9l / {k^ a (T- - T)) 1 / 2 , l 2 = 



[Hl/H^ a ']\ /2 g 2 /{n+a {T+ - T)) 1 / 2 , T± = T c ± 7 /a , and /i(i 4 ) is defined by 



(2u - ft - k)e- (u -^ 2 u- {1+ ^- l ^' 2 du = 0. (9) 
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Insight into the values of the coefficients gi and 52 can be obtained from the microscopic model originally studied 
by Ambegaokar et. al. for superfluid 3 He Q and extended to anisotropic superconductors by Sigrist and Ueda pl| 
and also by Samokhin The microscopic weak coupling analysis of [jl3| and |L4| indicates that for a specularly 
reflecting surface with normal n gi = if PjjA(fc) = A(fc) where PftA(k) = A(k — 2ft ■ kft) and that gi — 00 Jl4f] 
or gi « ceQT c <^Q,i ||l3|| (£o,i is the zero temperature coherence length) if P^jA(fc) = — A(fc). This is in agreement with 
experiments for conventional isotropic superconductors for which PfiA(k) = A(fe) for all surface normals with the 
consequence that g = 0. For a surface normal along a the above considerations imply that g 2 = and g\ — 00 or 
gi w a T c -Co,i |J] with Co,i = (^V/aoT-) 1 / 2 @. If we choose g x = c^T" Co,i Eqs. | and | indicate that H^ a * is 
equal to H® 2 to within 1/2 a percent for the temperature range we consider, an even stronger suppression of surface 
superconductivity will occur for any g\ > aoT c ~^o,i- Consequently gi > aoT c _ Co,i is equivalent in practice to g\ = 00. 
With these values for the coefficients gi two distinct H C3 curves are possible in UPt3, depending on whether 7 and 
5i + 64 are both positive or both negative. If 7 > and Si + 64 > then the bulk transition is to the 772 state (i.e. 
if = (0, 1)) for T* < T < T+ and to the 771 state for < T < T* . Since the 771 state is suppressed at the surface, 
surface superconductivity can only occur in the 772 state and the resulting temperature dependence of H^ a is shown 
in Fig. a (recall that for g 2 = 0, [H*f] 2 = 1.695[ff£j 2 )- On the other hand, for 7 < and 5x + S 4 < 0, the bulk 
transition is to the 771 state for T* < T < T c + and to the 772 state for < T < T* . In this case the behavior in Fig. |]b 
arises. 

The above discussion is relevant to specular reflecting surfaces which in fact does not account for the experimental 
results on UPt 3 (see below) . If diffusive scattering is considered then the 772 component will also be suppressed at the 
surface. This is taken into account by introducing g 2 7^ 0. In Fig. ||we have plotted H^ a /H% 2 versus temperature 
for gi > a T~( 0A and g 2 = a T+Co,2/10 where Co, 2 = («a /a T+) 1 / 2 . 

Keller has given a graph of AT(H) = T H a.a* (H)—Th<* 2 (H) versus H [[l0| that illustrates the temperature dependence 

of H^ a /H® 2 . It was observed that AT(H) initially increased with increasing field until it reached a maximum (at 
H « 0.25 Tesla) after which it decreased with increasing field. It can be seen that for specular reflection AT(H) is 
always increasing with H (see Fig. fa). The experimental behavior is best described by 7 > 0, Si + 64 > 0, and the 
diffusive scattering values aoT+( 0i2 /12 < g 2 < aoT+ £0,2/8 and gi > a T)TCo,i (as in Fig. 2). 

We have considered the simplest case in which the antiferromagnetic moment is pinned orthogonal to the surface. 
Other possible relative orientations of the surface and the antiferromagnetic moment are possible. Differing domain 
distributions at different surfaces may account for the different values of H C3 observed at these surfaces (note that 
different values of g 2 for different surfaces also can account for the observed H C3 values). A detailed study can be 
undertaken once a better understanding of the antiferromagnetic domain structure near the surface is achieved. We 
have also assumed that the whisker used by Keller et. al. can be described by the same model that describes annealed 
UPt3 samples. We note that Keller gives the H° 2 curve for this whisker p0[ . There is a strong curvature in H" 2 
at small fields and as the field increase above H « 0.25 Tesla (note this is approximately the same field at which 
AT(H) attains its maximum), H® 2 increases linearly with temperature. We expect that our model is qualitatively 
correct for this whisker. However, experiments on UPt3 samples with a well defined kink in H® will be required for 
a quantitative comparison. 

The model presented above also implies informative results for H°' a . For a surface normal along the c axis symmetry 
implies gi = g 2 in Eq. ^. It is therefore expected that H^ a will either exhibit a kink at T fa T* or be suppressed. 
Using the microscopic arguments presented earlier it is expected that for the E 2u and E\ g representations H^ a w H® 2 
while for the E\ u and E 2g representations H%' a > ff" 2 . 

A similar analysis for different models of UPt3 results in qualitatively different behavior than that of the E-KEP 
models. In the ID-REP model jl8| the order parameter is a three dimensional order parameter that transforms as a 
vector under spin rotations and as a one dimensional representation under spacial transformations. In this model the 
antiferromagnetic moment breaks the spin degeneracy of the order parameter. Since the order parameter transforms 
as a one dimensional representation under spacial transformations the surface free energy density takes the form <7|t7| 2 
for all surface normals (note that this model is exactly solvable for H^ a using Eqs. || and ^|). This surface free energy 
does not allow any anisotropy to occur between the different order parameter components. It is therefore difficult to 
reconcile this model with the experimental observations of Keller et.al. P,^0[. 

In the accidentally nearly degenerate AB model [[l9| the order parameter has two components of the same parity 
that transform as different one dimensional representations of the hexagonal point group (one as an A and the other 
as a B representation). The surface free energy for all surface normals takes the form £7^4 1 -0^4 1 2 + SbIV'bI 2 - This model 
is also exactly solvable for H C3 . Analysis indicates that the results for H^ a are similar to that of the E-KEP models 
with A = 0. However, a qualitative difference may occur in H^ 3 a between these two models. For H^' 3 a in the even 
parity AB models, microscopic analysis indicates that one component will be suppressed and the other will not be. 
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Consequently the resulting behavior of H^ a will correspond to that of one of the two curves in Fig. g. 

In conclusion, we have, for the first time, used experimental results on surface superconductivity to yield information 
about the symmetry of the superconducting order parameter. In particular the experiments of Keller et. al. |)|,|l0| on 
UPt3 can be understood if one component of the order parameter is strongly suppressed while the other component 
is weakly suppressed for a surface normal along the a direction. This anisotropy places a strong constraint on phe- 
nomenological models of superconductivity in UPt3. Further experimental investigations of surface superconductivity 
can place additional constraints on these models. 
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FIG. 1. H^ a (upper curves) and H£ 2 (lower curves) as a function of temperature for different model parameters for specular 
reflecting surfaces. The parameters have been chosen to agree with the experimental phase diagram for H£ 2 ■ a) 7 > 0, 81 +84 > 
b) 7 < 0, Si + 84 < 0. 

FIG. 2. The ratio JI°'° /Hc 2 f° r diffusive scattering at the surface. The rightmost curve corresponds to 7 > 0, Si + 84 > 
and the other curve to 7 < 0, 81 + 84 < 0. 
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